Cluster growth in far-from-equilibrium particle models with diffusion, detachment, 

reattachment and deposition 
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Monolayer cluster growth in far-from-equilibrium systems is investigated by applying simulation 
and analytic techniques to minimal hard core particle (exclusion) models. The first model (I), for 
post-deposition coarsening dynamics, contains mechanisms of diffusion, attachment, and slow acti- 
vated detachment (at rate e <C 1) of particles on a line. Simulation shows three successive regimes of 
cluster growth: fast attachment of isolated particles; detachment allowing further (et) ' coarsening 
of average cluster size; and t ' approach to a saturation size going like e _1 . Model II generalizes 
the first one in having an additional mechanism of particle deposition into cluster gaps, suppressed 
for the smallest gaps. This model exhibits early rapid filling, leading to slowing deposition due 
to the increasing scarcity of deposition sites, and then continued power law ((et) ) cluster size 
coarsening through the redistribution allowed by slow detachment. The basic (et) domain growth 
laws and e~ ' 2 saturation in model I are explained by a simple scaling picture involving the time 
for a particle to detach and diffuse to the next cluster. A second, fuller approach is presented which 
employs a mapping of cluster configurations to a column picture and an approximate factorization 
of the cluster configuration probability within the resulting master equation. This allows, through 
the steady state solution of the corresponding equation for a cluster probability generating function, 
quantitative results for the saturation of model I in excellent agreement with the simulation results. 
For model II, it provides a one-variable scaling function solution for the coarsening probability dis- 
tribution, and in particular quantitative agreement with the cluster length scaling and its amplitude. 



PACS numbers: 05.40.-a, 05.50.+q, 68.43. Jk 
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I. INTRODUCTION 

This paper is concerned with domain growth in far 
from equilibrium systems. This is a subject of increas- 
ing interest both for the wide range of behaviors and for 
the large number of applications, which range from phase 
separation in mixtures to island formation and coarsen- 
ing during deposition of a thin film or submonolayer 1,2 , 
among other systems. 

Our aim is to discuss a series of one-dimensional ex- 
clusion models with particle diffusion, reversible or ir- 
reversible attachment to clusters and deposition mecha- 
nisms that represent, for example, volume reduction ef- 
fects after cluster coalescence. Diffusion processes tend 
to bring these systems to equilibrium steady states, but 
pressure and other external influences may drive the sys- 
tem to new steady states. Though not usually related 
to a specific real problem, these one-dimensional models 
may reveal interesting features that help to understand 
more complex and realistic surface and bulk systems 3 ' , 
with the advantage of being more tractable both analyti- 
cally and numerically. We will discuss a series of plausible 
physical situations in systems with diffusion and mecha- 
nisms that drive them out of equilibrium, in order to un- 
derstand the details of domain growth and convergence 
to steady states, if it occurs. 

In the first model, hereafter called model I, a fixed frac- 



tion p of a one-dimensional lattice is randomly filled with 
hard core particles. The diffusion rates are r = 1 when 
they are free, i. e. when they have two empty nearest- 
neighbor sites, and r — e ~ e~ E ' T (where E is the re- 
lated energy barrier) when they have one occupied near- 
est neighbor site (Fig. la). For e < 1 (T -> 0), the aver- 
age aggregates' length grows as t 1 / 3 in a long time range, 
and eventually approaches saturation at ~ e -1 ' 2 with a 
slow i -1 / 2 decay (Sec. II). In the limit e — ► 0, this model 
is equivalent to the Ising model with Kawasaki dynamics 
previously studied by Cornell et al 5 , who focused on its 
zero temperature features. However, the dynamic rules 
are mainly motivated by the Clarke- Vvedensky model for 
thin films or submonolayer growth 6 , excluding the depo- 
sition processes. In the simplest versions of that model, 
an isolated adatom has to overcome an energy barrier E s 
to diffuse, while when it is attached to n nearest neigh- 
bors the energy barrier increases to E s + nE\,, where E\, 
is a bonding energy. This model and related ones were 
already intensively studied in two dimensions during the 
deposition process 6-8 , but a few works have considered 
the post-deposition coarsening dynamics 9 . 

Subsequently we will generalize the previous hard core 
dynamical system by introducing deposition of particles 
(see e. g. Refs. 11 and 12 ), but allowing deposition only at 
(empty) sites with at least one empty nearest neighbor 
(Fig. lb). In this model (referred to as model II), the 



domain coalescence, which generates larger vacancies be- 
tween aggregates, is followed by a density increase. The 
exclusion of deposition at single holes between clusters 
represent the geometrical frustration of real systems. In 
Sec. Ill, we will show that this model exhibits a t 1 / 2 
domain growth. This is among other results we have 
obtained by simulation studies, which are presented for 
both models I and II in Sees. II and III, respectively. 
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FIG. 1. (a) Diffusion and detachment processes of model I, 
with the corresponding rates, (b) Added deposition processes 
of model II, with the corresponding rates. The deposition at 
sites with two occupied neighbors (last process) is forbidden. 



Our models share some aspects with diffusion lim- 
ited coalescence models 15-17 and with fragmentation- 
aggregation models 13,14 . They are similarly described 
in terms of cluster or interval probabilities, and like the 
fragmentation models they are amenable to analytic in- 
vestigation based on an independent cluster approxima- 
tion (the independent interval approximation to the joint 
cluster length probability occurring in the master equa- 
tion) . We use this approach to explain properties of mod- 
els I and II, including distributions of cluster size (Sees. 
II and III). Further, a simple scaling picture can be devel- 
oped in order to describe the basic domain growth laws; 
we use this at the beginning of Sec. II. 



II. MODEL I: DIFFUSION, DETACHMENT AND 
REATTACHMENT OF PARTICLES 

A. Processes 

The model studied in this section has the particle hop- 
ping processes depicted in Fig. la. Isolated particles hop 
symmetrically on a chain at unit rate ("diffusion"), while 
a single particle with a left hand / right hand neighbor 
can hop to an empty right / left neighbor with rate e (de- 
tachment). So clusters evolve by detachment and reat- 
tachment of particles. The model is of exclusion type: no 
site can accommodate more than one particle. 

This model is clearly particle-conserving, so density p 
is fixed. The case e ^C 1 is of particular interest since, as 
reported in the simulation studies below and explained 
in the following subsection, very large clusters emerge. 



B. Simulations 

We simulated model I in one-dimensional lattices of 
length L = 8000. This length is sufficiently large to en- 
sure that finite size effects are negligible, as shown by 
comparisons of some results with data from lattices with 
L = 16000 (particularly for the smallest values of e this 
comparison is essential). 

Initially, the lattice is randomly filled with a density 
of particles p. We simulated three values of the density, 
p = 0.1, p = 0.5 and p = 0.9, which are representative 
of the range of intermediate densities, i. e. densities not 
too small {p « 0) nor too large (p sw 1). For p = 0.1 and 
p = 0.9, we considered several values of the diffusion rate 
e ranging from e = 10 _1 to e = 10~ 3 , and for p = 0.5 we 
performed simulations until e = 10~ 5 . 

The sequence of characteristic behaviors of model I, as 
shown by simulation results, are: (i) early fast attach- 
ment of isolated particles to each other to form clusters; 
(ii) an intermediate regime in which detachment sets in, 
allowing further coarsening; (iii) finally, there is a dif- 
fusive approach to a saturated state where the clusters 
have a large steady mean size that depends on e. 

The three regimes are well separated at small e. This is 
illustrated in the plot of log 10 d versus time log 10 t, shown 
in Fig. 2, where d is the mean size of clusters of two 
or more particles; d is given in terms of the probability 
P t (m) that an arbitrarily chosen cluster has size (or mass) 
m at time t by 



£m=2 mP t M 
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FIG. 2. Typical time evolution of the average cluster length 
d in model I, with three different regimes. Data in the plot 
were obtained for p = 0.5 and e = 10~ 4 . 



The early time dependence of d in region (i) (at small 
e) starts with a characteristic increase with rate propor- 
tional to a high power of e, and then crosses over to a 
form allowing data collapse in terms of the reduced time 
variable et, as shown in Fig. 3. 
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In region (ii), d increases as 



d~t a . 



(2) 



The apparent exponent a e ff , defined as the local slope of 
the logd x logi plot, was calculated numerically. a e // is 
shown in Fig. 4 as a function of (et)~ for three different 
values of e. It appears to approach the value a = 1/3 
in the limit of small e and correspondingly large t, which 
is consistent with the prediction of a simple scaling de- 
scription (next subsection). 

Fig. 5 shows the diffusive approach of the mean clus- 
ter size to its saturation value doo . This approach is well 
described by d — d^ — -Spi-, for t — ► oo, with C constant. 
The dependence on e of the saturation value d^ is illus- 
trated in Fig. 6 for p = 0.5. The least squares fit in Fig. 
6 gives 



0.72e- 1/2 + 1.93, 



(3) 



in which the dominant (proportional to e -1 ' 2 ) and the 
sub-dominant (additive constant) terms were estimated. 
Like Eq. (2) with a = 1/3, this result follows from the 
analytic work in Sec. II C. 



FIG. 3. Average cluster length as a function of the scaling 
variable log 10 (et), for model I with p = 0.5 and e = 10~ 3 
(squares), e = 10 -4 (triangles) and e = 5 x 10~ 5 (crosses). 
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FIG. 5. Average cluster size at long times t as a function 
of l/i 1/2 , for p = 0.5 and e = 10" 3 in model I. The dashed 
line is a least squares fit of the data for 10 < t < 10 . 



FIG. 4. Effective exponents a, defined as local slopes of the 
logd x log* plots, for model I with p = 0.5 and e = 5x 10~ 5 
(below, blue), e = 2 x 10 -5 (medium, green) and e = 10 _J 
(up, red). 
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be replaced by an effective rate e = eP m i g , where P m %g is 
the (migration) probability that a freely diffusing particle 
does not return to the detachment site before diffusing 
the distance ~ d to the next cluster. In other words (see 
Fig. 7b), this is the probability that a free particle at 
position y = 1 does not return to the origin (y = 0) 
before a time of order d 2 , which is the typical time for 
diffusion along a distance d. Considering that 



Q(y,t) 



y 



(4ttL>) 



1/2 



-3/2 



exp 



ADtJ 



(5) 



is the probability that the first passage of a random 
walker at point y occurs at time t 18 (D is the diffusion 
coefficient), P m i a is given by 



P„ 



<p 



Q(l,t)dt~d~ 



(6) 



FIG. 6. Saturation value of the average cluster length as a 
function of e~ ' 2 in model I. The dotted line is a least squares 
fit of the data. 
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C. Theory 

The characteristic results just described have been in- 
terpreted by simple heuristic scaling arguments and by 
detailed analytic studies starting from the master equa- 
tion and employing an independent interval approxima- 
tion. This second approach is capable of yielding cluster 
length distributions and their evolution. 

To begin with, we focus on the asymptotic cluster size 
at small e. This asymptotics occurs in the regime where 
lone particles are rare, and those that are present are in 
the process of reattaching themselves to the clusters they 
came from or to a neighboring one. The second case pro- 
vides the sharing which sets the mean cluster size d. At 
densities of order p ~ 1/2, the cluster size is roughly of 
the order of cluster separation (see Fig. 7a). Thus the 
equilibrium of detachment time (1/e) and time of diffu- 
sion to a neighboring cluster (~ d 2 ) gives the observed 
saturation result 



,-i/2 



(4) 



This argument can be generalized rather obviously to 
explain the t~ 1 ' 2 approach to saturation. 

A more interesting application is the explanation of 
the early cluster size growth law (Eq. 2). Here, unlike 
the saturation just described, the cluster separations are 
such that the detached particle is likely to return and 
reattach many times before it eventually diffuses to the 
next cluster (Figs. 7a and 7b). Its likelihood of returning 
to the origin means that the detachment rate e needs to 
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FIG. 7. (a) Detachment of a particle (in grey) at the border 
of a cluster, with rate e. In model I, the mean cluster length 
is d and the typical cluster separation is, asymptotically, also 
of order d for densities not too small nor too large, (b) A free 
particle (in grey) immediately after its detachment from the 
right cluster. 



In terms of the effective rate e, the required time for 
a particle to transfer to the next cluster is of order 1/e. 
Thus, the time required for doubling the size d of a cluster 
by successive gain and loss of particles is d 2 /e ~ d 3 /e. So 
the cluster growth proceeds according to 



—d d 

dt ~WN) 



and hence 



(et) 



1/3 



(7) 



(8) 



This explains the behavior seen in the simulations (Figs. 
3 and 4). The situation is analogous to domain scaling in 
Ising chains where, with Kawasaki dynamics 5 , spins split 
off from domain edges and migrate across to increase the 
domain size by one lattice unit. 

We turn next to the more powerful analysis starting 
from a version of the master equation, which can provide 



a full description of the process. This is more easily set 
up by reformulating the process using a column picture, 
in which a column of height m represents a cluster of 
size to, and then the original detachment and diffusion 
processes correspond to those shown in Fig. 8. Since one 
cluster has two edges but corresponds to a single column, 
the one-particle detachment rate in the column picture 
is 



7 = 2e. 



(9) 
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FIG. 8. (a) Example of particle-hole configuration on a line 
and the map (dashed arrows) into a column problem, (b) The 
processes of particle detachment from clusters, with rate 7, 
and of free particles diffusion, with rate 1, in the correspond- 
ing column picture. 



We denote by Pt(m) the probability that a randomly 
chosen cluster (equivalently, column) has size to at time 
t. Then the gain/loss from in/out processes provides the 
following master equation, in an independent interval ap- 
proximation in which joint probabilities are factorized: 



P t+1 (m) - P t (m) = A r . 



Pt(rn- l)0(m-l) 



7 ^2 Pt(m') + 1 x P t (m' = 1) 



m'>2 



P t {m + 1) [ 7 0(ro + 1 - 2) + 1 x 6 m ,o] 

j P t (TO)[ 7 0(TO-2) + lx<5 mjl + 

7 J2 Pt(m') + lxP t (m' = l)}. 



(10) 



The corresponding equation for the generating function 

G t ( s )^]Tp t (TOK" (ii) 



G t +i{s) = G t (a) l + sa{t) + ±-j-a(t) 

l s 

+s ( 7 -l)P t (l) + 

[ 7 P t (0) + (l- 7 )P t (l)]-Ip t (0), (12) 

where a(t) = 7 [1 - P t (0)] + (1 - 7)P*(1). ft is easy to 
check probability and mass conservation using Gt (0) and 

G'M- 

The steady state distribution P{m) and generating 
function G(s) resulting from Eq. (12) are given by 



G(s) = (7 - say 1 [ 7 P(0) - s(l - 7 )P(1)] 



P(m) 



P(1),TO> 1, 



and 



P(l) = AP(0), 



(13) 



(14) 



(15) 



with A = 7 [1 - P(0)] [1 - (1 - 7 )P(0)] . So the steady 
state cluster size distribution is exponential. The mean 
size of multi-particle clusters (Eq. 1) and the mean mass 
(to) = X)m=o m P ( m ) are t nen obtainable in terms of 
P(0), as is the density p. So, in particular the mean 
cluster length d = (to) can be found in terms of p. The 
result simplifies at small 7 (small e) to 



d = 1 - 1 ' 2 



-1/2 



(1-P) 



1/2 , (p/(l- P ) + 'i) 



p 1 1/2 , ( P /(l-p) + 3) 



2(1-P) 



(16) 



is 



where the dominant and the first sub-dominant terms are 
shown. This form is consistent with the scaling result (4) 
and is in very good agreement with the simulation result 
(3), including the sub-dominant constant term: Eq. (16) 
gives d = 0.7071 . . .e" 1 / 2 + 2 for p = 0.5 (see also Fig. 
6). In the same limit, this is also, apart from a numerical 
factor, the characteristic size in the exponential cluster 
mass distribution. 



III. MODEL II: DIFFUSION AND DEPOSITION 
OF PARTICLES 

A. Processes 

Model If is a generalization of model I, different only 
by having the deposition processes depicted in Fig. lb, 
in addition to the diffusion and detachment processes of 
Fig. la. This makes the model non-particle-conserving, 
which leads to continued coarsening and other scaling 
properties and crossover. 



B. Simulations 

The characteristic behavior of model II, as exhibited 
by simulation results, is as follows. For initial densities 
not too near p = 1, there is: (i) an early regime of rapid 
hlling, due to deposition, and cluster evolution due to 
both processes; (ii) an intermediate regime where depo- 
sition slows because of the scarcity of deposition sites 
due to the increased density - the exclusion constraint of 
course applies. The slow detachment process allows re- 
distribution of particles, opening up new deposition sites 
and allowing the continually slowing coarsening (with no 
saturation as p < 1). 

Fig. 9 shows simulation results for the evolution of the 
mean cluster size d. That plot shows that d(t) is well 
fitted by the form 



d(t) = Bt 1 ' 2 (l + cr 1 ' 2 



(17) 



In Fig. 10 we show the ratio between the estimates of 
the amplitude B and e 1//2 for several values of e. Those 
results give 



BQ-fe 1 /", 



(18) 



with negligible corrections to scaling, where b = 0.252 ± 
0.002. This result is in accord with theoretical analysis 
given in the next subsection, including the estimate of 
the amplitude b. 
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FIG. 9. Long time scaling of the mean cluster length (d) in 
model II. 



1 /9 

of the scaling variable (ei) ' and, at very long times, it 
converges to 1 as 



i- P ~(e*r 1/2 - 



(19) 



This form is also in agreement with the theoretical anal- 
ysis below. 
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FIG. 10. Ratio between the estimated amplitude B of the 
scaling of average cluster length (Eq. 17) and e 1 , as a func- 
tion of log 10 (1/e), in model II. 
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The dependence of the evolving density on time t and 
rate e has also been studied. The simulation results 
shown in Fig. 11 imply that the density is a function 



FIG. 11. Scaling plot of the particle density p in model II, 
for e = 1CF 2 (squares), e = 5 x 10~ 3 (triangles) and e = 1CP 3 
(crosses) . 



C. Theory 

The characteristics presented in the foregoing subsec- 
tion can be interpreted using an analytic investigation 
along the lines of the detailed discussion given in Sec. 
IIC. 

We have to include the effects of the extra deposition 
process, which leads to the decrease of the total number 
of clusters and of the number of holes between the clus- 
ters as time increases. On the other hand, the length L 
of the line in which particles are deposited and diffuse is 
kept constant. Consequently, in order to adopt the col- 
umn picture of Sec. II. 3 (see Fig. 8), it is necessary to 
consider that the length Lq of the corresponding column 
problem decreases in time (these lengths are related as 
Lq = L — M, where M is the total mass or total number 
of particles, for periodic boundary conditions). 

The evolution equation here is written for cluster num- 
bers as 



N(m, t + 1)- N(m, t) = L (A m + B m ) 



(20) 



where the diffusion contribution A m is given in Eq. (10) 
and the deposition contribution is 



B m = P t (0)[29 (to - 2) P t (to - 1) - 
20 (to - 1) P t (to) + S mtl P t (0) - 2<5 m , ] 



(21) 



The length of the lattice in which the column problem is 
defined varies due to deposition as 



L (t + 1) - L (t) 

L (t) 



= -P t (0)[2-P t (0)] 



In these equations, the cluster probability is 

N(m,t) 



Pt(m) 



L (t) 



(22) 



(23) 



They preserve conservation of probability, but mass is no 
longer conserved. 

The resulting equation is similar to ones occurring in 
coalescence models 15-17 . From this we expect that a 
large time and small e limit discussed subsequently is 
equivalent to the model in Ref. 15 . Our approach, which 
exploits the generating function method, becomes equiv- 
alent, in the scaling limit, though in a conjugate space, to 
continuum approximations used in the coalescence stud- 
ies of Refs. 15 and 17 . 

Now the generating function (Eq. 11) satisfies 

G t+1 (s)[l - Pt(0)} 2 - G t (s) = 

(s 1) \G t (s) (a(t) - 1) + ( 7 - l)Pt(l) + 2p t (0)l + 

L V s / si 

P t (0) [2(s - l)G t ( S ) - sP t (0) + 2 (P t (0) - 1)] . (24) 

In the right hand site of Eq. (24), the first term cor- 
responds to diffusion processes and the second one to 
deposition processes. 



Because deposition slowly fills the system, we expect 
the configurations to coarsen and presumably to go into 
some scaling asymptotics where mass scales with some 
power of t, and Pt(m) and Gt(s) each become one- 
variable scaling functions. So we look for a long time 
scaling solution of the above equation. 

At long times, the finite difference Gt+i(s) — Gt(s) in 
Eq. (24) can be taken as a derivative. The scaling vari- 
able will be some combination of t (large) and u = 1 — s 
(small), the latter because large cluster sizes arise from 
structure in G t (s) at s w 1. The variable u is actually 
conjugate to to (see below). Coarsening will correspond 
to the scale of to as t z , with some power z, in which case 
the one- variable form will be 



G t (s) = u a f(ut z ). 



(25) 



with some function /. Normalization requires a = 
and /(0) = 1. In the scaling limit, the relationship of the 
generating function to the probability P t (to) requires the 
latter to be of the form 



with 



1 /TO 



f(x) - / g(y)e- xy dy. 
Jo 



(26) 



(27) 



It turns out that the consistent scaling solution has 
g(0) = 0, so the \/t z contribution to P*(0) vanishes, leav- 
ing a leading term of lower-than-scaling order, 



*M = % 



(28) 



where c is a constant. Eq. (24) leads to the dynamical 
exponent 



z = 1/2 



(29) 



and to the following equation for the one- variable scaling 
function: 

xf'(x) - 2cf(x) - 2 1 x 2 f{x) + 2c = 0. (30) 

Even without solving Eq. (30) we can infer that 

c/2 



P(0) 



t ' 



and 



(to), ~ 7 1/2 * 1/2 ~ e 1/2 t 1/2 



l-p t ^ 7 - 1 / 2 t- 1 / 2 ^e- 1 /V 1 / 2 . 



(31) 



(32) 



(33) 



These hold in the long time scaling limit we have intro- 
duced and agree with the observed simulation results in 
Eqs. (17) , (18) and (19). 



Eq. (30) can be formally solved for the scaling function 
f(x) by using the variable £ = x 2 and considering the 
function f(x)x~ 2c . The result is 



m 



1 Jo 



dv. 



(34) 



The large x expansion of f(x) and Eq. (27) provides the 
small y expansion of g(y): 



9(y) = ^(^J, 



(35) 



where 



Gin) = t < C + ll^ + m \ -i r ^. (36) 



m— 



This confirms that g(0) = 0. The cluster distribution has 
the following form, in terms of the odd function Q: 



Pt(m) 



1 M "' 



m*(t) \m*(t)J' 



(37) 



where to* = ("ft) ~ (et) . It explains the scaling 

1 /2 

variable (et) used to collapse simulation data in Fig. 
11. 

The conditions that Q must be non-negative and nor- 
malisable are satisfied with c = 1/2 in Eq. (28), which 
leads to 



G(u) = \e~^ 



(38) 



The mean cluster mass (cluster length in the original 
problem) is easily obtained as 



V^(7*) 



1/2 



(39) 



Considering relation (9), the amplitude of cluster length 
scaling is B = V2tt£ 1/2 « 2.5076 1 / 2 . It quantita- 
tively agrees with the result obtained in simulations (Sec. 

iiib"). 



IV. CONCLUSION 

We studied two one-dimensional exclusion models with 
particle diffusion, reversible or irreversible attachment to 
clusters and deposition mechanisms. 

In model I, starting from a randomly filled lattice, only 
particle diffusion is allowed, with small detachment rates 
e for particles at the edges of the clusters. Simulation 
results show an initial regime with formation of small 
clusters, a regime of cluster size growth as d <~ t 1 / 3 and a 
regime of cluster size saturation at d ~ (T 1 ' 2 . These re- 
sults can be explained using heuristic scaling arguments. 
The analytical treatment of the master equation with an 
independent cluster approximation for joint probabilities 



distributions predicts a saturation cluster size in quanti- 
tative agreement with numerical data. 

Model II generalizes model I in having also particle 
deposition: this is allowed only at empty sites with one 
or two empty nearest neighbors. Simulation results show 
continuous coarsening with a t 1 ' 2 increase of the aver- 
age cluster size and an increase of the density with t -1 ' 2 
corrections. These scaling forms are justified by ana- 
lytical investigations again using an independent cluster 
approximation, which provides good quantitative agree- 
ment with the simulations. 

We expect that the models presented above and the 
combination of different methods to explain their scal- 
ing behaviors can be used to understand further non- 
equilibrium systems. Of particular interest would be 
the extension of theoretical methods (e. g. scaling ap- 
proaches) to two-dimensional systems such as adatom is- 
lands on surfaces or the extension of the one-dimensional 
models to include other mechanisms that drive the sys- 
tems to new non-equilibrium steady states or which lead 
to anomalous coarsening. 
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